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Abstract 

We compute the effective action for stationary and spatially constant magnetic fields, when 
coupled anomalously to charge neutral fermions, by integrating out the fermions. From this the 
grand canonical partition function and potential of the fermions and fields are computed. This 
also takes care of magnetic field dependent vacuum corrections to the grand canonical potential. 
Possible applications to neutron stars are indicated. 
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I. INTRODUCTION 


The study of magnetized nuclear matter systems finds application in physical systems 
at various length scales and temperatures: It has been proposed that magnetic fields of the 
order of 10^® G can be induced in non-central heavy-ion collisions [l|. Magnetic fields of 
similar stre^th can also be found in the interior of highly magnetized neutron stars, called 
magnet ars [2|, yj. 


Neutron stars are extremely compact objects with densities similar to, and greater than, 
that of nuclear matter. These are stable systems and assumed to be in thermal equilibrium 
at zero temperature, since the thermal energy is much less than the Fermi energy of the 
particles in the star’s interior. Matter with similar densities is also formed in relativistic 
heavy-ion collisions. Although this type of matter is stable for brief periods, due to the 
much higher temperatures, it is also assumed to be in thermal equilibrium 


Under the assumption of thermal equilibrium, the appropriate description of these sys¬ 
tems is a statistical description based on the grand canonical ensemble with the grand 
canonical partition function and closely related potential as the central quantities. Quite 
generally, therefore, the computation of the grand canonical potential for nuclear matter at 
hnite temperature and density in the presence of strong magnetic helds is highly relevant. 
The most convenient formalism for this is finite temperature feld theory as described in 
many textbooks. Here we follow the formalism as presented in 


his quantity using 
^ hnite temperature 


Due to its high relevance there were a number of recent studies of_t 
a variety of techniques and focussing on a variety of contributions. In 
effects for baryons in strong magnetic helds were studied by computing the ehective action 
for the electromagnetic held. More recently explored the thermodynamical properties 
of strongly magnetised fermions with emphasis on the issue of pressure anisotropy induced 
by the magnetic held. In both these works the anomalous coupling and possible vacuum 
contributions to the grand canonical potential were not taken into account i.e., the “no sea” 
approximation was made. It is important to emphasise that the anomalous coupling referred 
to here is not just the radiative corrections for charged elementary particles, which one ex- 
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pects to have a minor effect [8|, but it also arises in effective theories that treat structured 
baryons with an associated magnetic dipole moment as point particles. To compensate for 
this it is required to include an anomalous coupling to the dipole moment of such particles in 
their description. This is, of course, the case for charged and charge neutral baryons, but it 
should be clear that it is of particular importance for charge neutral ones, such as neutrons, 
since it is their only interaction with the electromagnetic field. The role of this coupling 


for magnetised fermions was first explored in [9j and has since been used extensively to 
investigate magnetised systems. Indeed, in a recent study the role of an anomalous 
coupling in determining the bulk properties of a gas of charged and charge neutral fermions 
was explored quite extensively, again emphasising the role of pressure anisotropy. However, 
this study also did not include any possible vacuum contributions to the grand canonical 
potential. 


To appreciate the crucial role that the thermodynamic potential plays, recall that at equi¬ 
librium the thermodynamic potential is the quantity that is minimised and from which all 
thermodynamic quantities are computed. It should therefore be clear that for a full under¬ 
standing of the equilibrium properties of these systems it is imperative to have as complete 
computation as possible of the grand canonical partition function. In particular one should 
recognise that the vacuum energy may also yield magnetic field dependent contributions. 


An equivalent way of phrasing this, in the language of quantum field theory, is as fol¬ 
lows: To find the full magnetic field dependence of the grand canonical potential, one has to 
compute the effective action of the electromagnetic field by integrating out the fermions in 
the presence of a chemical potential. This automatically takes care of all the contributions, 
including vacuum, of the fermions. When doing this, one treats the electromagnetic field as 
classical, i.e., one does not integrate over it. The computation of the relevant determinants, 
even with classical electromagnetic field configurations, is in general intractable. However, 
these determinants can be computed when restricting to stationary and spatially constant 
field configurations, which is the general mean field approximation made in all the works 
mentioned above. This was also the approach first taken in jo], but not including any 
anomalous couplings. Of course, when following this approach one has to face ultraviolet 
divergences, which we address in some detail here. 
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As far as we can establish, no compntation of the effective action for the electromag¬ 
netic held in the presence of anomalous couplings has been done. The aim of this paper is to 
address this, particularly in the case of neutrally charged fermions that couple anomalously 
to the electromagnetic held. The motivation for studying this particular case is that it is the 
most relevant scenario for neutron rich nuclear matter in strong magnetic helds as found in 
neutron stars. Furthermore, the extension of these results to charged, anomalously coupler 
fermions is straightforward. Apart from these motivations, current literature 7|, [ill, [12 
still contain several controversies, especially around the issue of anisotropic pressure. We 
indicate a possible resolution of this controversy here, but, since this is not the main aim of 
this article, a more careful and detailed analysis will be carried out elsewhere. 


Once the ehective action has been derived, one can extract the full magnetic held de¬ 
pendence of the grand canonical potential, which can then in turn be used to study the 
equilibrium and, at zero temperature, the ground state properties of the system. 


II. FORMALISM 


A system of magnetised neutrons is described by 




( 1 ) 


where 'ifj is the neutron held operator and | [ 7 ^, 7 ^"] are the generators of the Lorentz 

group. In this paper we will use the Heaviside-Lorentz units and the convention of Ref. 13 1 
for the 7 -matrices (for exact expressions of quantities and their units in this conventions 


please see Ref. 


nl) 


As already mentioned, the coupling between the dipole moment of the baryons and the 
electromagnetic held tensor is reminiscent of the coupling arising from the anomalous cor¬ 
rection to the electron magnetic dipole moment. This is a higher order quantum correction 
due to radiative corrections to the photon-electron vertex. Since neutrons are not funda¬ 
mental particles this coupling does not represent such a quantum correction, but rather an 
ehective coupling between the neutron’s magnetic dipole moment arising from its internal 
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structure and the magnetic field. However, in keeping with convention we will still refer to 
this coupling as the “anomalous” magnetic moment. The coupling strength gb is therefore 
the neutron’s magnetic dipole moment in units of fermi. 


Our aim is calculate the grand canonical partition function and grand canonical poten¬ 
tial, or more simply the grand potential, of this system using the Feynman path integral 
^proach, as described in Ref. The grand canonical partition function is expressed as 

3 


Z = j [dvr] / [d(j)] exp 




dr I d^x — gAf{(j),7T) 


L^o 


( 2 ) 


with /3 = {kBT)~^ the inverse temperature times Boltzmann’s constant, r = —it the an¬ 
alytic continuation of time in the complex plane, Ti the Hamiltonian density, and J\f the 
conserved density. The functional integration over 0 represents all fields, while that over n 
their conjugates. 


We assume that the magnetic field is a stationary and spatially constant background 
field pointing in the 2 ;-direction. Based on these assumptions we choose the gauge field 
to = (0,0,i?x,0). This choice simplifies —in ([T]) to and reduces the 

anomalous coupling to 

- = gb^ ■ B, (3) 

where S = cr 0 I 2 , cr are the Pauli matrices, and I 2 the 2 x 2 identity matrix. Since the 
magnetic field is considered to be a classical background field only the fermion fields are 
integrated over in Z. 


The Hamiltonian density is 


( 4 ) 


'ip - C = 'ip [-ij ■ V - gbZ; ■ B + m\'ip + 


( 4 ) 


and the conserved density in ([2]) is the neutron number density. Thereby Z reduces to 


Z = I [id'ip^] [d'lp] exp / dr I cPx ip ( 7 '^^—h ^7 ■ V -|- gbZ ■ B — m — ) tp 


d 


B^ 


.(5) 


ar ' ” ' ■ ' r-, j -r 2 

From here on we will momentarily ignore the contribution of the Maxwell term since this 
only adds a constant factor. Since the addition of the anomalous term does not affect the 
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commutation relations of the baryon fields, the method of Ref. j5| can be used to evaluate 
the functional integrals. Thereby Z reduces to the determinant of the matrix D which does 
not contain any held operators, 


so that 




D ^ -i/3 [(-iw„ + m) - ■ V - ■ B + m)] , 


( 6 ) 


In Z = In Det D = Tr In D = \n [I3‘^[{ujn + inf+ , (7) 

\ n k 

where Un = (2?7, + 1) T are the Matsubara frequencies. The u is the spectrum of the 
Hamiltonian (jl]) which is known from Ref. [l^ (and references therein): 


A) 



+ ^2) 


( 8 ) 


where 


A = ±1 distinguishes the different orientations of the neutron dipole moment, and 


• k]_ = + ky, perpendicular to B = Bz, given the choice of 

The sum over n is computed in and thereafter ([7]) reduces to 

(Pk 


\^z = vY^ f 

\ ^ 


R2 


, , ^ [^co + In (1 + + In (1 + - /3R , 

(ztt)'* 2 


( 9 ) 


where we have again included the Maxwell contribution of the magnetic held. The grand 
potential is dehned as 

\nZ 


f2(/i,T) = 


and its various thermodynamic relations are 


vp ’ 


( 10 ) 


s 


dn\ 


N 


V V 


dn\ 


( 11 ) 


Due to the presence of the magnetic held the spherical symmetry of the ground state is 
broken and the pressure is generally anisotropic. The longitudinal P\\ and transverse P± 
components of the pressure are given by 

dil 


p = -n, p± = -n + B 


dB' 


( 12 ) 
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( 13 ) 


In the limit of /3 ^ cx), so that T —)■ 0, for // 7^ 0 (and fi > m) fl becomes 

^ j -uj + {uJ- ©[a; - fj]} + —. 

Using the following relation between the various thermodynamic quantities, dE = TdS — 
PdV + fiN, the energy density of the system is 

^=y = J2j + + (14) 

Therefore the hrst term of the integrand in flT^ corresponds to the zero-point energy of the 
vacuum, while the second term is the contribution from all the hlled positive energy particle 
states (since T = 0 the chemical potential /i is equal to the neutron’s Fermi energy). 


Since we are interested in the dependence of hi on the magnetic held B, all terms that 
are independent of B are irrelevant. In particular this implies that the only relevant quan¬ 
tity is the difference between the grand potential for B 0 (magnetised) and B = 0 
(unmagnetised), which we expand in a power series in B. The difference in the contribu¬ 
tions from the occupied positive energy neutron states is hnite and calculable for a hxed 
number of neutrons and we will refer to it as A(/i,i?). The difference in hi for magnetised 
and unmagnetised neutrons is therefore written as 

00 

Q{B) - Q{0) = A{fi, B) + J2 «2n(A)S'", (15) 

n=0 


where A is some (large) momentum cut-off that in the end will be taken to inhnity. The 
expansion coefficients a 2 n are determined by expanding the single particle energy (|H]) for 
5^0 


J2‘^{b) = j2^{o) + 


siKB^ 


(5A)2-4(P + m2))S^ 


+ 0(5®). (16) 


The hrst term in the expansion is just the zero-point energy for unmagnetized neutrons and 
therefore will be cancelled in flT^ . In the remaining terms the k± integrals can be performed 
without any divergences appearing, while a momentum cut-oh A is introduced for the k^ 
integrals. Dropping terms of order and higher yields. 


= / (0, 


(1 - log[4]) 
47r^ Stt^ 


+ 


2 2 
9b^ 

47J-2 


log 


m 

A 


B^ + 


f/6log[2] 

247r2 


+ 


ffb 

247r2 


log 


m 

,A. 


(17) 

5^ 
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Furthermore B < 1 and therefore higher order terms in B become irrelevant. We 

note the appearance of quadratic and logarithmic divergences when A —)■ oo. As we do 
not have any physical renormalisation point, we remove the inhnities through a minimal 


subtraction scheme 


15[ |. In addition, we must keep in mind that as here represents the 


full grand potential of matter and helds, it must yield, upon setting gi, = 0, the contribution 
of the Maxwell term ^B^ and free fermions, i.e., A(/x, S) = 0, which implies that 02 = ^ 
and a 2 n = 0, Vn 7 ^ 1 at this point. Then the magnetic held dependent part of the grand 
potential for magnetised neutrons is 


^ / _ glm^ (1 - log[4]) \ _ glB'^ log[2] 

2 \ / 247r2 

In the zero temperature limit of /? —)■ cxo flTSl) becomes 

B^ f g^m^ {1 - \og[4]) g^B^ \og[2]\ 


(18) 


(19) 


where p = y and 

^ = j ■ (20) 

A 

Not surprisingly we note that the fermion contribution to the effective action for the elec¬ 
tromagnetic held modihes the coefficient in front of the B^ term, which simply implies that 
the vacuum is polarised due to the presence of the matter helds. 


III. DISCUSSION 

The zero-point energy is the inhnite contribution from the hlled negative energy (antipar¬ 
ticle) states that form the vacuum. For most dense matter systems this “sea” of antiparticles 
is ignored since only the relative changes in the hlled positive energy states are assumed 
to have physical consequences. Generally this is a good approximation for dense nuclear 
matter system since the high value of the Fermi energy, up S> —m, blocks any low energy 
vacuum excitations due to the hlled positive energy states. Figured] demonstrates the ehect 
of the vacuum contribution, using the free value of ^ 1.913p,N = 90 ~ 0.0305 fm, to 
the thermodynamic potential and conhrms its relative smallness. 











B [lO'^G] 


FIG. 1. (Color online) The contribution of the zero-point energy to the last term in (|19p . 
Correspondingly it should be compared to 1, which represents the contributions from the magnetic 
field. 


It should, however, be kept in mind that since we treated the electromagnetic held as 
classical, photon loop corrections to gt,, and therefore also medium effects, have not been 
accounted for. These may signihcantly alter the value of gb and therefore also the effect 
of the vacuum contribution. Indeed, there is good experimental evidence that the neutron 
dipole moment may increase with density: in Co ppe r isotopes the dipole moment changed 
by about 50% over an atomic mass range of 10 1^. This suggests that the contribution 
of the vacuum energy may be much more relevant if all radiative corrections, and thus 
medium effects, are more carefully accounted for, which is one of the main motivations 
for computing these corrections here. One can, of course, simply treat gb as a parameter 
and investigate its potential impact on the equilibrium properties of these systems. This 
suggests considerable quantitative and ^alitative effects and may even, if gb is made large 
enough, lead to a ferromagnetic phase 1^. These are, however, currently speculations that 
need to be scrutinised carefully in future. 


A hnal observation we make here relates to the rather controversial issue of the anisotropic 
pressure. Recall that at equilibrium —O = —Of + with Of the fermion contribution. 
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is minimised or, alternatively, the entropy is maximised at fixed temperature and density. 
From this one can compute the equilibrium value of the magnetic field as 


B=^ = M, ( 21 ) 

dB ' ^ ^ 

where M is the magnetisation of the fermion gas. Substituting this in flT^ shows that at 
equilibrium 

P\\ = P± = + — • ( 22 ) 

This suggests that an anisotropic pressure may be a signal that the system is out of equi¬ 
librium and, since B ^ M, the presence of free currents and a dynamo effect. These are 
merely hints that warrant further investigation elsewhere. 

IV. CONCLUSION 

In this article we have computed the full magnetic field dependence of the grand potential 
for magnetised neutron matter coupled anomalously to the magnetic field. This includes 
vacuum effects but not yet radiative photon corrections. We did this by computing the 
effective action of the electromagnetic field, treated as a classical field, by integrating out the 
fermions. The result can be used to compute the properties of these systems at equilibrium 
and this may shed light on the controversial issues of anisotropic pressure and the possibility 
of a ferromagnetic phase transition. 
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